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ABSTRACT 

We present an investigation of quasar colour-redshift parameter space in order to 
search for radio-quiet red quasars and to test the ability of a variant of the KX quasar 
selection method to detect quasars over a full range of colour without bias. This is 
achieved by combining IRIS2 imaging with the complete Fornax Cluster Spectroscopic 
Survey to probe parameter space unavailable to other surveys. We construct a new 
sample of 69 quasars with measured bj — K colours. We show that the colour distri- 
bution of these quasars is significantly different from that of the Large Bright Quasar 
Survey's quasars at a 99.9% confidence level. We find 11 of our sample of 69 quasars 
have signifcantly red colours (6j — K ^ 3.5) and from this, we estimate the red quasar 
fraction of the K < 18.4 quasar population to be 31%, and robustly constrain it to 
be at least 22%. We show that the KX method variant used here is more effective 
than the UVX selection method, and has less colour bias than optical colour-colour 
selection methods. 

Key words: (galaxies:) quasars: general, techniques: photometric, galaxies: active, 
infrared: galaxies, X-rays:diffuse background 



1 INTRODUCTION 



IWebster et al.l (|l995l ) compared the B — K colours 
of Parkes Half-Jansky Flat-Spectrum Sample (PHJFS) 
iDrinkwater et all [l993) q uasars to those of Large Bright 
Quasar Survey (LBQS) (|Hewett. Foltz. fc Chaffed Il995l ) 
quasars. The PHJFS quasars were found to have a signif- 
icantly broader and redder distribution of B — K colours 
than LBQS quasars. This disparity should not exist accord- 
ing to th e unified model o f Active Galac t ic Nu clei (AGN) 
galaxies (lAntonuccilll993h . IWebster et al.l |l99^ concluded 
that there exists a populatio n of red quasars that at that 
time had not been discovered. IWebster et al.l |l99^ investi- 
gated this possibility using a simple model of LBQS selec- 
tion effects, and demonstrated that the LBQS was biased 
against the detection of red quasars. The bias was caused 
by the blue magnitude limit of the LBQS, because in a 
blue-flux limited survey, red quasars need to be intrinsically 
brigh ter than blue quasars to be detected (jWebster et al.l 
The bias causes fewer red quasar detections, because 
the luminosity function of qua sars is steep at the bright end 
iBovle. Jones, fc Shank'3ll99ll ). Other surveys with a blue- 
flux limit are similarly biassed, independent of the quasar 
selection method. Using their simple model, IWebster et al.l 
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(|l995h estimated that the red quasars missed by previous 
surveys constitute as much as 80% of all quasars. 

Red quasars are defined by a colour criterion. The best 
colour to use for this criterion is that made from filters sep- 
arated by the largest wavelength range possible, because 
anomalously redder quasar colours stand out more. For the 
datasets used in this paper, the colour covering the largest 
wavelength range is U — K. We used the next largest colour 
bj — K, as far more objects, ~ 3 times as many, have a 
measured bj — K colour compared to the number with a 
measured U — K colour. The additional advantage of having 
used bj — K is that it is possible to compare the colours of the 
quasar sample presented here with the B — colour used in 
[Glikman et al (200j) and the g — K colour of lHopk ins et al.l 
120041). Throughout this paper we define a r ed quasar by the 
criterion, bj - K ^ 3.5 (|Gregg et al.ll2002l ). because when 
describing the optical-infrared region of quasar spectra with 
a power-law, ^(A) cx A", this colour corresponds to an effec- 
tive optical- infrared spectral index of a ^ 1. 

The existence of red quasars, other than in radio- 
selected quasar samples, over most of the bj — K colour 
range of the PHJFS quasars, satisfying t he above cri- 
terion , has already been d emonstrated |Hopkins et al.l 
20041: [Richards et al. I l2003l : iGlikman et al.l |2004| ). The 



Glikman et al. 1 l|2004l ) red quasar sample contains the red- 
dest quasars found, with B — K colours of 5 ^ B 
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K ^ 8. iHopkins et al] (|2004l ): iRichards etU] (|2003l ) 

found red quasars in the Sloan Digital Sky Survey 
(;SDSS) (|York et al.ll200ol ). and iGlikman et all (lioM) found 
them in a combination of the 2-Micron All Sky Survey 
(2MASS) jSkrutskie et al.l |2006| ) and an Automatic Plate 
Measuring CATalogue (APMCAT) of POSS/U KST sky sur- 
vey plates (jlrwin. Maddox. fc McMahon|[l994h . 

Despite confirming the existence of red quasars, current 
observations have not explained the mechanism or mecha- 
nisms responsible for their existence. Dust obscuration by 
dust at the quasar redsh ift, in either the hos t galaxy or the 
quasar accretion disk ( Webster et al.l Il995l l , was the first 
mechanism proposed for quasar reddening. IWebster et al.l 
l| 19951 ) proposed this mechanism as a correlation between 
B — K colour and redshift was not observed. The result of 
reddening a composite optically selected quasar spectral en- 
ergy distribution (SED) with dust at the quasar redshift, 
was found to be consistent with the dust mechanism, be- 
cause t he reddening resulte d in the SED of a red PHJFS 
quasar (|Webster et al.lll995l). 

As an alternative, iBenn et al.l | l99H) sugg ested that 
the B — K spread in the iWebster et al.l (jl995f ) sample is 
partly or entirely caused by host galaxy starlight. Quasars 
with fl at-radio-spectra a re typically hosted by later galaxy 
types l|Benn et al.l 1 19981 ). and the old stellar population 
of late type galaxies is a stronger em itter in K than 
B, reddening quasars (jBenn et al.l Il998l ) . An analysis in 
iMasci. Webster, fc Francis! (|l998l ) of the host galaxy con- 
tri bution to the r e d B — K colours of the PHJFS quasars 
m IWebster et all l|l995h was inconclusive, but suggested 
host galaxy contribution was insufficie nt to cause the red 
B — K colours observed. Further work in lMaddox fc HewettI 
l|2006t ) confirmed that host galaxy starlight can definitely 
redden quasars, particularly low luminosity and low redshift 
resolv ed quasars. In accordance with Maddox fc Hewett 
l|2006l ). in this paper we use the criteria that luminous un- 
resolved quasars, which are not located at low redshift, are 
unlikely to be redd ened by host galaxy c o ntribu tion. 

Alternatively, ISerieant fc Rawhngd (|l99d ) proposed 
that a synchrotron component of quasar emission cre- 
ated a K-band excess, reddening quasars. An anal- 
ysis of 157 PHJFS quasars and 12 LBQS quasars 
found the shape of the SEDs was consistent with both 
line-of-sight dust and synchrotron emi ssion mechanisms 
Francis. Whiting, fc Webster 2 000 ). In iBarkhouse fc Halll 



iMalhotra. Rhoads. fc Turner! ll 19971 ) observed a similar phe- 



20011 ) it was suggested that because any synchrotron com- 



ponent to quasar emission is weaker in radio-quiet quasars 
than in radio-loud quasars, any synchrotron emission will 
only contribute to a small fraction of an y reddening. Based 
on the quasar sample contained therein, IBarkhouse fc Hall! 
l|200ll ) suggested that any radio-quiet quasar redder than 
bj — K ^ 3.7 is too red to be the result of synchrotron 
emission, and is most likely the result of dust obscuration. 
Throughout this paper we use the criterion that any red 
quasar that is also radio-quiet, is unlikely to be reddened by 

synchrotron emission. 

Finally, IMa lhotra. Rhoads. fc Turneil l|l997h suggested 
that a small subset of red quasars was caused by the lens- 
ing of a quasar by a dusty galaxy. Not only does the dust 
in the lens redden the quasar but the dusty galaxy po- 
tentially magnifies the host galaxy starlight contribution 
l|Gregg et al.!l2002h . This last mechanism was proposed after 



nomenon to IWebster et al.l (1993); with radio-selected lenses 
having redder colours than optically selected lenses. 

Whichever mechanism creates red quasars, they are pos- 
sibly an evolutionary stage. After obser ving low-ion i sation 
Broad Absorption Line (BAL) quasars, iLacv et all (|2002D 
suggested a model of quasar evolution where a quasar occu- 
pies a thick dusty torus after forming, eventually the torus 
dissipates leaving behind an optically blue quasar. Alter- 
nately, quasars may emerge from dusty star burst galaxies, 
residing in the dust stirred up by mergers jUrrutia et al.l 
[2005). Either evolutionary sequence would naturally cre- 
ate a situation where quasars are reddened by dust at the 
quasar redshift. A third evolutionary path has been pro- 
posed where red quasars are an evolutionary link between 
Ultra-Luminous Infrared Galaxies (ULIRGs) and optically 
selected quasars (^Chen et al. 2006). Determining the frac- 
tion of red quasars and any dependence on redshift, lumi- 
nosity, radio fiux, etc. will indicate if red quasars are an 
evolutionary stage, and which evolutionary path red quasars 
belong to. 

Depending on the mechanism, red quasars 
might explain the quasar contribution to the X- 
Ray Background (XRB) . Using a theoretical model, 
iMadau. Ghisellini. fc Fabian! l|l994 ) showed that a ratio of 
two to three times as many dust obscured to unobscured ac- 
tive galactic nuclei (AGN) galaxies can adequately account 
for the existence of the XRB. The model also accounts for 
the observed spectrum and the source counts in the soft 
and hard X-ray bands. The latest theoretical models of 
the XRB r equire a smaller ratio, lying somew here between 
0.6 to 1.5 (|Gilh. Comastri. fc Hasinge j [2OO7I ) . If it can be 
shown that red quasars are the result of dust obscuration, 
and that the red quasar fraction is within 37.5% to 60%, 
then red quasars are the most probable source of the quasar 
contribution to the XRB. 

The existence of red quasars has various implications, 
the main one is that it can no longer be assumed that ex- 
isting quasar samples are representative of the entire quasar 
population. The statistics of a sample refiects the proper- 
ties of the sub-set of a population from which the sample 
is drawn, previous quasar samples, such as the LBQS, are 
biassed to the inclusion of blue quasars; therefore, the statis- 
tics of such samples predominantly refiects the blue quasar 
population. Use of existing quasar datasets requiring a sam- 
ple that is representative of the entire quasar population is 
therefore subject to review. Direct examples of this are that 
the luminosity function of quasars and the change in number 
density with redshift need to be re-evaluated while including 
red quasars. 

In order to learn more about red quasars, we need to 
spectroscopically observe quasars over their entire colour 
range. To achieve this, future surveys will have to select 
quasar candidates for spectroscopic follow-up using a tech- 
nique that is not biased against either red or blue quasars. 
One candidate for a suitable selection technique is th e KX 
method, developed in lWarren. Hewett. fc Folt j (|2Q00'). The 
KX method utilises the power-law nature of quasar SEDs 
at long wavelengths combined with morphological classi- 
fication. Stars experience a turnover in their SED in H- 
band while quasars follow a power-law, allowing discrimi- 
nation between stars and quasars. The stellar morphology 



Finding red quasars and testing the KX method 3 



of quasars discriminates between quasars and galaxies. Ob- 
servationally, this requires constructing a colour-colour plot 
of all objects with a stellar morphology, using an optical 
colour and a colour that straddles the H-band. 

In this paper we explore regions of the quasar colour- 
redshift space unavailable in previous work, and test a KX 
method variant. We aim to find the boundaries of the quasar 
colour-redshift parameter space that is occupied. In doing 
so, we will improve the estimate of the red quasar fraction 
and possibly shed light on the mechanism that creates red 
quasars. In testing a KX method variant, we intend to verify 
if it is a viable method of selecting quasars, without bias, 
from the entire quasar colour-redshift parameter space. We 
intend to apply a KX method variant to our quasar sample, 
and then contrast it with established methods of selecting 
potential quasars for spectroscopic follow-up. 

To achieve these goals we use a sub-sample of the 
quasars ident ified in the Fornax Clus ter Spectroscopic Sur- 
vey (FCSS) ijPrinkwater et al]l200Cl l. a spectroscopic sur- 
vey of all extended objects to hj = 19.8 and point sources 
to bj = 21.5. The sub-sample consists of all quasars that 
match to an object in Ks-band imaging to a depth of K = 
18.4, obtained using the Infrared Imager and Spectrograph 
2 (IRIS2) instrument on the 3.92-m Anglo- Australian Tele- 
scope (A AT). The Ks-band imaging is deeper than 2MASS 
allowing us t o explore parts of the bj — K colour space 
unavailable in iGlikman et all l|2004l ). which was limited to 
fe./ - ^ 5 by the K = 14.5 magnitude limit of 2MASS. 
Furthermore, the quasar sample use d covers a larg e r rang e 
of redshift than was examined in iHopkins et al. (|2004h : 
iRichards etHI ()2003h . which were limited to z ^ 2.2. Be- 
cause the FCSS used no selection criteria, our quasar sam- 
ple is bias free, and ideal for testing the ability of the KX 
method to select quasars from a K magnitude limited survey. 
Additionally, most of the FCSS objects have a U magnitude 
as well as a b,j magnitude, allowing the KX method to be 
directly contrasted with traditional methods for optically 
selecting quasars. 

The structure of this paper is as follows. In Section [2] 
the datasets used in this paper are described. Section |3] de- 
tails the detection of red quasars in the quasar sample and 
an analysis of the quasar sample. This analysis involves a 
comparison of the bj — K colours of the sample quasars 
to LBQS quasars and a determination of the effect of our 
selection on the quasar sample constructed here. Section [3] 
concludes with the calculation of the fraction of red quasars. 
The effectiveness of the KX method is assessed in Section |4] 
and following this are our conclusions in Section [5] 



2 PHOTOMETRIC AND SPECTROSCOPIC 
DATA 

The sample of quasars used here and their photo- 
metric information were obtained by combining three 
datasets. Astrometry and optical photometry were ob- 
tained from the Automatic Plate Measuring CATa- 
logue (APMCAT) of POSS/UKS T sky survey plates 
l|lrwin. Maddox. fc McMahonI Il994h . IRIS2 imaging was 
used to obtain the K-band photometry used here. The nec- 
essary spectroscopic identifications were provided by the 
Fornax Cluster Spectroscopic Survey (FCSS) catalogue of 



iDrinkwater et al.l ||200C| ). The datasets were combined by in- 
dependently matching the astrometry of K-band catalogue 
objects and FCSS objects to that of APMCAT catalogue 
objects. In the rest of this section, we describe the datasets 
in more detail. 

The FCSS catalogue is a blind spectroscopic survey 
of Fornax objects in the APMCAT catalogue. The FCSS 
is an ideal source of spectroscopic identifications because 
the lack of target selection precludes biasing. In the FCSS, 
spectroscopy was carried out on all extended sources to 
bj ^ 19.8 and all point sources to bj ^ 21.5. Not every 
source was spectroscopically observed, and the fraction of 
observed sources is referred to as the spectroscopic com- 
pleteness. The FCSS obtained a redshift and spectroscopic 
identification for more than 90% of the objects observed, re- 
ferred to as the redshift completeness. In most cases where 
a spectroscopic identification could not be obtained, neither 
was a redshift because of the poor spectrum quality. Be- 
cause the redshift completeness was consistently high for 
the entire FCSS, we have combined the redshift complete- 
ness with the spectroscopic completeness, and throughout 
the paper the term 'completeness' refers to this combina- 
tion. The completeness of the FCSS is 96% for bj < 20.5, 
82% for 20.5 ^ bj <^ 21 and 36% for 21 < fej ^ 21.5 
sources. We account for the completeness by appropriately 
weighting quasars during our analysis in later sections. 

IRIS2 was used in December 2001 to obtain Ks-band 
imaging using 60s exposures. The raw data were processed 
using the Observatory Reduction and Acquisition Control 
project - Data Reducti o n pipeline (ORAC-DR ) outli ned in 
lEconomou et al.l (|l999l ): Ijenness fc Economoul (|l999l 'l. The 
reduced IR IS2 imaging dataset w as then analysed with 
SExtractor (|Bertin fc Arnouts|[l996l ') to obtain a photomet- 
ric catalogue of corrected isophotal magnitudes. Corrected 
isophotal magnitudes were used instead of PSF or small 
aperture magntiudes, because PSF and small aperture mag- 
nitudes minimise any possible extended contribution to the 
flux of sources in the IRIS2 imaging. The Ks-band cata- 
logue magnitude zero-point was then calibrated, by match- 
ing point sources in the Ks-band catalogue to the 2MASS 
point source catalogue. Calibrating also facilitated the con- 
version from Ks magnitudes to the K-band magnitudes used 
throughout this paper. After calibration, the magnitude 
limit of the K-band photometry was found to be 18.4 mag- 
nitudes. 

In the FCSS, quasars were spectroscopically identified 
as objects with broad permitted lines with a measured fu ll 
width half-maximum > 1100 km s"^ (|Mever et al.l 120011 ). 
After matching the IRIS2 imaging catalogues and the FCSS 
quasar identifications to the APMCAT object positions us- 
ing blind matching, to a radius of 10"; a catalogue of all the 
necessary photometry was created, and all quasars within 
both the IRIS2 imaging K magnitude limit, and the FCSS 
spectroscopic identifications bj magnitude limit were identi- 
fied. This resulted in a sample of 69 spectroscopically iden- 
tified quasars with U, bj and K magnitudes. All 69 of these 
quasars are APMCAT point sources; therefore, our quasar 
sample is limited to K ^ 18.4 and b,j ^ 21.5. Of these 69 
quasars, 62 have an R magnitude as required by our KX 
method variant. The missing R magnitudes are caused by 
the different magnitude limits in the bj and R bands of the 
APMCAT catalogue, 22.5 and 21 magnitudes respectively. 
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Figure 1. A plot of bj K colour as a function of absolute bj 
magnitude, with (right) and without (left) uncertainties. Quasars 
with a red, bj — K ^ 3.5, colour arc denoted by circles and blue 
quasars with solid triangles. The bj — K colour uncertainty is the 
result of combining in quadrature the photometric uncertainties 
in the bj and K ma g with the quasar variability bj magnitude 
structure function in lHook et al. 

3 RED QUASARS 

Combining the APMCAT, IRIS2 imaging and FCSS 
datasets resulted in a sample of 69 spectroscopically identi- 
fied quasars. In Section [3. II we examine the bj — K colours 
of the quasars to determine which quasars are red. We in- 
vestigate the uncertainties in the data to ascertain if the 
uncertainties caused blue quasars to appear red. Section[3]2] 
compares the colours of the constructed quasar sample to 
the colours of a sub-set of LBQS quasars to characterise the 
'reddness' of the quasar sample. Then in Section [3.41 we ex- 
amine selection effects and estimate the red quasar fraction. 
Finally, Section[3]3]examines the relative bj—k colours of our 
quasar sample, and compares them to the relative colours of 
the LBQS sub-set. 

3.1 Detection of red quasars 

We show the bj — K colours of our quasar sample in Figure 
[T] 11 quasars are red because they satisfy bj — K ^ 3.5, 
corresponding to ~ 0.16 of the 69 quasars with bj and K 
magnitudes. These red quasars have colours, 3.5 ^ bj — K ^ 
5.6, w hich occupy the colour space between I Glikman et al.l 
l|2004) red quasars and LBQS, SDSS red quasars. Using 
the All-Sky Optical Catalo gue of Radio/X-Ray Sources in 
iFlesch fc Hardcastld (|2004l ) . we verified the existence (or 
lack thereof) of a radio-counterpart for all the red quasars 
discovered here, to a flux limit of ImJy. We found that 5 of 
these quasars had no radio-counterpart and the remaining 
6 were not in the catalogue, indicating they too are radio- 
quiet. We concluded that all of the red quasars in F igure [1] 
are radio-quiet quasars. In lMaddox fc HewettI l|2006l ) it was 



Figure 2. bj - K colour vs redshift, z, for the 69 quasars presented 
here. Blue quasars (bj — k < 3.5) are solid triangles and red 
quasars {bj — K ^ 3.5) are circles. Using redshift bins of size z 
= 0.7 (reference image in bottom right corner) centered at points 
0.07 apart from z = 0.3 to z = 2.8, the number of quasars and 
median bj — K colour is plotted. The median bj — K colour is 
plotted as a dashed line and the number of number of quasars in 
each bin, with a separate axis on the right hand side, is plotted 
as a dotted line. 

demonstrated that the host galaxy contribution to quasar 
flux is maximal in the K-band for resolved, low-luminosity 
quasars ai z < 1. bj — K is plotted against redshift in Fig- 
ure [2] and a dotted line tracks the number of quasars as 
a function of redshift, peaking at z « 1.5. The dashed line 
in Figure [2] traces the median quasar colour with redshift, 
demonstrating a dependence of the observed bj — K colour 
on redshift. Section |3]3] examines the dependence of the ob- 
served bj — K colour on redshift. As the red quasars found in 
Figure[T]are all luminous unresolved APMCAT sources, and 
the majority are luminous quasars with 2; > 1 in Figure [2] 
any host galaxy contribution to the bj~K colours of the red 
quasars will be minimal. If a significant host galaxy contri- 
bution was present in the red quasars identified in Figure [TJ 
then the bj — K colours of the red quasars would become in- 
creasingly redder as both luminosity and redshift decreased. 
In Figures [1] and [2] a clear trend of increasingly red bj — K 
colours with decreasing luminosity and redshift is not appar- 
ent. The 11 red radio-quiet quasars detected here are most 
likely the result of reddening by dust located at the quasar 
redshift. 

We investigated the uncertainties in the observed bj — K 
colours, to determine if the uncertainties caused blue quasars 
to appear red. Uncertainty in the bj — K colours in Fig- 
ure [T] has two sources, photometric uncertainties in the 
b,j and K magnitudes and the effect of quasar variability 
on multi-epoch observations. The photo metric uncertainty 
in the bj magnitudes was obtained from iDrinkwater et al] 
(|2000l ) where it was found to be 0.1 mag. For the K-band 
magnitudes the uncertainty was calculated by SExtractor 
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Table 1. Completeness limits and associated weighting applied 
to a quasar as a function of bj magnitude. 



bj magnitude 


completeness 


weight applied to quasar 


bj < 20.5 


96% 


1.04 


20.5 < bj ^ 21 


82% 


1.22 


21 < bj ^ 21.5 


36% 


2.78 
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Figure 3. Results of scattering the bj — K colours of 51 blue 
(fej — K < 3.5) quasars with Gaussians scaled by the photometric 
uncertainties (dotted line), quasar variability structure function 
(short dashed line) and both combined in quadrature (solid line) . 
Plotted are histograms of the percentage of 10,000 iterations that 
resulted in a given fraction of the blue quasars being scattered 
so as to have red (fej — K ^ 3.5) colours. The long dashed line 
marks the fraction of red quasars found here. 

during the data reduction. To account for the efTect of quasar 
variabihty we estimated the uncertainty in the bj photom- 
etry, due to it being observed at an earlier epo ch to the K 
band p hotometry, using the structure function of lHook et al.l 

\Amtj \ = [0.155 + 0.023(AfB + 25.7)]AfSi,i^(l + z) (1) 

where Atobs is the time difference between the multi-epoch 
observations in the observed frame. In doing this we effec- 
tively move the bj photometry to the epoch of the K pho- 
tometry. To evaluate the structure function we used Atobs = 
25 years and calculated t he absolute bj mag nitudes using 
k{z) = -1.25/og(l + (|Mever et all 1200 ll ). = 0.3, 
r^A = 0.7 and Ho = 72 km s"^ Mpc"^ The bj ~ K colour 
uncertainties in Figure [1] are the combination of the pho- 
tometric uncertainties and the structure function result in 
quadrature. Based on Figure [TJ some of the red quasars 
found here are potentially blue quasars, that appeared red 
because of uncertainties in their bj — K colours. 

We explicitly tested if all the red quasars observed here 
are the result of photometric uncertainties and/or quasar 



variability by scattering the measured bj — K colours of the 
51 blue quasars {b,j — K < 3.5) by a Gaussian, scaled so 
that 1-sigma equals the bj — K uncertainty, and calculat- 
ing the resultant fraction of red quasars; after weighting the 
quasars by the completeness of the FCSS, using the weights 
in Table [T] The solid line in Figure [3] shows the result of 
10,000 attempts at scattering the bj — K colours of the blue 
quasars; a histogram of the percentage of the 10,000 itera- 
tions that resulted in a given red quasar fraction. The red 
quasar fraction observed here, after weighting the quasars 
using the weights in Table[TJ is ~0.22. Based on the percent- 
age of 10,000 iterations that resulted in an equal or greater 
red quasar fraction, we concluded from Figure |3] that we 
could reject at the 99.9% confidence level that our observed 
red quasar fraction is the result of photometric uncertainties 
and/or quasar variability. 

The solid line in Figure |3] has a peak at a non-zero red 
quasar fraction, suggesting that the bj — K uncertainties 
are reddening our quasar sample. To investigate the cause 
of this potential reddening we again scattered the bj — K 
colours of the 51 blue quasars, this time scaling the Gaus- 
sians by the bj — K photometric uncertainties and the effect 
of quasar variability separately, and the result of 10,000 it- 
erations are shown in Figure [3] as a dotted and dashed line 
respectively. From Figure |3] we concluded that quasar vari- 
ability was potentially reddening our quasar sample as the 
dashed line in Figure |3] peaks at a non-zero red quasar frac- 
tion; and that the photometric uncertainties did not redden 
our quasar sample. 

3.2 Comparison to LBQS quasars 

The sample of 69 quasars was compared to a sample of 538 
LBQS quasars that matched objects in the 2MASS survey. 
Histograms of both sets of quasar b,j — K colours were plot- 
ted in Figure m where the bj — K colours for our 69 quasars 
were weighted by the completeness of the FCSS, using the 
values shown in Table [1] The histograms in Figure |4] sug- 
gest that the quasars presented here have a different bj — K 
distribution to the LBQS quasars, covering a larger range 
and extending noticeably redder. For the weighted bj — K 
colours, the average, standard deviation, excess kurtosis and 
skewness are bj — K = 2.9, a = 0.9, —0.48 and 0.46 mag- 
nitudes respectively. The 538 LBQS quasars have an av- 
erage, standard deviation, excess kurtosis and skewness of 
bj ~ K — 2.8, a = 0.5, —2.7 and 0.04 magnitudes respec- 
tively. Plotting cumulative frequency distributions of the 69 
measured bj — K colours and the 538 LBQS quasars in Fig- 
ure (5] a Kuiper test showed the two distributions to be dif- 
ferent at the 99.9% confidence level. The statistics of the two 
b,j — K distributions are consistent with our sample having 
a redder, broader distribution than the LBQS quasars. 

The concern in making this comparison is that any dif- 
ference between the quasar sample constructed here and the 
538 LBQS quasars may be the result ofbj — K uncertainties. 
To investigate this we scattered the bj — K colours of the 
LBQS quasars by the uncertainties in their b,j — K colours, 
using Gaussians scaled such that l-cr equals the bj — K un- 
certainty, and tested if the two distributions were still dif- 
ferent using a Kuiper test. The LBQS bj — K uncertainties 
were determined by adding in quadrature the bj magnitude 
uncertainties in the LBQS, K magnitude uncertainties in 



6 R.J. Jurek, M.J. Drinkwater, P.J. Francis and K.A. Pimhhlet 




1 2 3 4 5 6 

bj - K (mag) 



Figure 4. Histograms of quasar fraction as a function of bj — K 
colour. The solid line is the weighted bj — K distribution of the 69 
quasars in our sample, weighted using the values in Table [l] The 
dashed histogram is the distribution for a sample of 538 LBQS 
quasars matched to 2MASS. The distribution of the 69 measured 
hj — K colours extends noticeably further and redder than the 
distribution of LBQS bj — K colours. 
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Figure 5. Cumulative distribution function of bj — K colours. 
The solid line is for the 69 quasars in our sample, weighted using 
the values in Table [T] The dashed line is for a sample of LBQS 
quasars matched to 2MASS. A Kuiper test of the two distribu- 
tions shows them to be different at the 99.9% confidence level. 



2MASS and the bj — K uncertainties due to quasar vari- 
ability. The bj — K uncertainties due to variability were ca l- 
culated usins the structure function of iHook et all l|l994h . 
outlined in Section lSTT] with Atobs = 11 years. Scattering the 
LBQS b,j — K colours 10,000 times, a Kuiper test showed 
the distributions of the scattered LBQS bj — K colours and 
our quasar samples bj — K colours were different 99.3%, 
71.7% and 24.9% of the time, at the 95% (2-cr), 99% and 
99.7% (3-(t) confidence levels respectively. Based on the re- 
sults of scattering the LBQS b,j — K colours, we concluded 
that bj — K uncertainties were unlikely to cause the ob- 
served difference between our quasar sample and the LBQS 
quasars. 

3.3 Relative bj — K colours 

As the observed bj — K colour is dependent upon redshift, 
there may be red quasars in our sample that in the observed 
frame failed to meet the bj — K ^ 3.5 colour criterion; in- 
stead being classified as blue quasars. We investigated this 
in Figure |6] by plotting relative bj ~ K colours, the quasar 
bj — K colours minus the median bj — K colour for the quasar 
redshift, essentially applying a first order K-correction. We 
measured the median bj — K colour, that was subtracted 
from both our quasar sample and the LBQS quasar sample, 
at a given redshift from the LBQS quasar sample; because, 
in a quasar population where red quasars occupy a tail that 
extends away from blue quasars, the median bj ~ K colour 
will be dominated by the bj — K colours of blue quasars. 
All the quasars in our sample, that are red based on their 
observed bj — K {bj ~ K ^ 3.5) colour have a relative colour 
such that bj ~ K 0.7, as such we define this as the relative 
bj — K colour criterion for a red quasar. One quasar that 
is blue, based on its observed bj — K colour, would have 
been classified red using its relative bj — K colour. This 
quasar demonstrates that observed frame colour criteria are 
potentially unsuitable for different iating betw e en bl ue and 
red quasars, as commented on in iHall et al] (|2006l ). It is 
beyond the capability of this paper to determine whether 
observed or relative bj — K colours are better suited to dif- 
ferentiating between blue and red quasars, because of the 
small size of our quasar sample. 

To investigate the effect of using relative colours on our 
comparison, we compared the relative b,j — K colours of our 
quasar sample to the relative colours of the LBQS quasars. 
We constructed the sample of relative LBQS colours by sub- 
tracting the median LBQS bj — K colour, at a given redshift, 
from each LBQS quasar. In Figure [7] we plotted the rela- 
tive bj — K distributions of both our sample and the LBQS 
quasars. In Figure[7]our samples relative bj — K distribution 
appears to be broader and redder than the LBQS relative 
b,j — K distribution; and, compared to Figure |4l the small 
amount of excess blue quasars in our sample compared to 
the LBQS decreased, while the excess of red quasars com- 
pared to the LBQS increased. Our samples relative b,j ~ K 
colours have an average, standard deviation, excess kurtosis 
and skewness of bj — K = 0.2 and a = 0.9, —0.6 and 0.5 
magnitudes respectively. The relative colours of the LBQS 
quasars have an average, standard deviation, excess kurtosis 
and skewness of bj - K = -0.01, a = 0.43, -2.8 and 0.01 
magnitudes respectively. Plotting cumulative frequency dis- 
tributions of the relative bj — K colours of the 69 quasars 
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Figure 6. Figure |2] re-plotted with the median colour for a given 
redshift subtratcted to give relative bj — K colours. Quasars with 
blue observed bj — K colours (bj — K < 3.5) are plotted as solid 
triangles and quasars with red observed bj — K colours {bj — 
K ^ 3.5) are plotted as circles. A dashed horizontal line marks 
the relative bj — K criterion proposed for red quasars, relative 
bj — K ^ 0.5. Four quasars with blue observed bj — K colours 
have red relative bj — K colours, demonstrating the dependence 
of red/blue quasar classification on redshift when using observed 
bj — K colours. 

in our sample, and the relative hj — K colours of the LBQS 
quasars in Figure [H a Kuiper test showed the two distribu- 
tions to be different at the 99.9% confidence level. The statis- 
tics of the two relative bj — K distributions are consistent 
with our quasar sample having a broader and redder relative 
b,j — K distribution than the LBQS quasars. Based on the 
differences between Figures |4] and [T] and our statistical anal- 
yses here and in Section [3. 2 1 we determined that whether we 
used observed or relative quasar colours, the results of com- 
paring our quasar sample to a sample of LBQS quasars is 
the same, that our sample has a broader, redder bj — K 
distribution than the LBQS quasars. The additional infor- 
mation gleaned from our comparison of the relative bj — K 
colours is that the small excess of blue quasars in our sam- 
ple, when compared to the LBQS, appears to be caused by 
K-correction effects, and is not actually an intrinsic prop- 
erty of the bj — K distribution. Because the excess of blue 
quasars in our quasar sample, when compared to the LBQS 
quasar sample, is probably caused by K-correction effects, 
the small excess of blue quasars is much less significant than 
the excess of red quasars when comparing our quasar sample 
to the LBQS sample. 



3.4 Limits on the red quasar fraction 

To understand the effects of our sample selection on the bj 
and K magnitude parameter space, we plotted bj and K 
magnitudes in Figure [9] with completeness limits and lines 




relative bj - K (mag) 



Figure 7. Histograms of quasar fraction as a function of relative 
bj — K colour. The solid line is the weighted distribution of our 
samples 69 quasars relative bj — K colours, weighted using the 
values in Table[T] The dashed histogram is the relative bj — K dis- 
tribution for a sample of 538 LBQS quasars matched to 2MASS. 
The distribution of our quasar samples relative bj — K colours 
extends noticeably further and redder than the distribution of 
LBQS relative bj ~ K colours. 
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Figure 8. Cumulative distribution function of relative bj — K 
colours. The solid line is for the relative bj — K colours of the 
69 quasars in our sample, weighted using the values in Table [T] 
The dashed line is for the relative bj — K colours of a sample 
of LBQS quasars matched to 2MASS. A Kuiper test of the two 
distributions shows them to be different at the 99.9% confidence 
level. 
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Figure 9. bj vs. K for the 69 quasars presented here. Blue quasars 
(bj - k < 3.5) are sohd triangles and red quasars (bj - K ^ 3.5) 
are circles. Lines of constant bj — K are plotted as dashed lines, 
indexed as they increase from the bottom right corner. Solid hor- 
izontal lines mark the regions of the FCSS completeness as a 
function of bj magnitude, and are labeled accordingly. 



Figure 10. Plot of the observed bj — K (left) and K (right) distri- 
butions of our quasar sample, and the distributions predicted by 
our model of quasar density in the bj — K vs. K parameter space. 
The solid lines are the observed distributions and the dashed lines 
are the distributions predicted by our model. The model is a rea- 
sonable match, in both shape and scale, to the observed bj — K 
and K distributions. 



of constant bj — K overlaid. The plot in Figure [9] revealed 
where in bj and K magnitude parameter space the number of 
quasar detections declines; and, by extension revealed where 
the bj — K distribution will drop off due to selection effects. 
The upper left region of Figure [5] is where we expect the 
least quasar detections due to the combination of complete- 
ness, and the magnitude limits of the bj and K photometry, 
referred to collectively as selection effects. A few quasars 
are found with bj — K > 4 and only one is found satisfying 
K < 17 and bj > 20.5. As the decrease in observed quasar 
detections coincides with where it is expected to decrease 
because of selection effects; we conclude from Figure [9] that 
the decrease in quasar detections, redder than bj — K = 4, 
is a result of sample selection and not an intrinsic property 
of our quasar sample. This is consistent with the fact that 
redder qu asars than thos e foun d here have been observed, 
such as in iGlikman etall (|2004l ). 

The analysis of our sample selection supports the hy- 
pothesis that the bj — K distribution extends further and 
redder than observed here, as such our observed red quasar 
fraction is a robust lower limit of the actual red quasar frac- 
tion for K ^ 18.4. Weighting the quasars by the complete- 
ness of the FCSS, using the weights in Table [1] we calcu- 
lated the observed red quasar fraction and found that a ro- 
bust lower limit for the actual red quasar fraction is --^0.22. 
To estimate the red quasar fraction for K ^ 18.4, we first 
assumed that the bj — K colour distribution is indepen- 
dent of K magnitude. Our second assumption in estimat- 
ing the red quasar fraction, was that the true extent of the 
bj — K distribution of radio-quiet quasars should match that 
of radio-loud quasars; therefore, we assumed that the true 
extent of the bj — K distribution for our sample is bj — K 



= 8, matching the extent of t he bj — K distr i butio n of the 
PHJFS radio-loud quasars in IWebster et al] |l99i). Using 
these two assumptions, we postulated a model for the num- 
ber density of quasars in bj — K vs. K parameter space, then 
used maximum likelihood estimation to find values for the 
model parameters; and finally integrated this model to esti- 
mate the red q u asar f racti on. As Figure]^ and t he results in 
[Richards etlo] l|2003t ) and iHopkins et al.l (|2004l ) are consis- 
tent with a.bj — K colour distribution, where red quasars oc- 
cupy a tail extending away from blue quasars, for the bj — K 
component of our model we used a Gamma distribution. For 
the K component of our model we have assumed a power- 
law, 10^''. Normalised by the number of quasars in our sam- 
ple, the resultant model for the number density of quasars 
in bj — K vs. K parameter space is, 

N{K,bj - K) ^ew'^^e-^^^'^ibj - Ky-^ (2) 

and maximum likelihood estimates of the model parame- 
ters are 6 = 46.65, /? = 0.39, a = 0.36 and 7 = 8.56. Due 
to limited data, we were unable to test the validity of our 
model in describing the data statisitically; therefore, we as- 
sessed the model validity by plotting the observed bj — K 
colour and K magnitude distribution with that predicted by 
our model in Figure 1101 and comparing them. Comparing 
the observed and predicted distributions in Figure 1101 we 
concluded that the model was suitable for estimating the 
red quasar fraction; because it reproduced the overall shape 
and scale of both the K and bj — K distributions. Integrat- 
ing our model over 15 ^ K ^ 18.4, for the bj — K ranges 
Q ^b,j — K ^ 8 and 3.5 ^ bj — K ^ 9,, we calculated a red 
quasar fraction of 31%. Using the observed b,j — K distri- 



Finding red quasars and testing the KX method 9 



bution of our sample, we therefore estimate the K ^ 18.4 
red quasar fraction to be 31%, and robustly constrain it to 
be greater than 22%. If the majority of red quasars are the 
result of dust reddening, our estimate of, and lower bound 
on , the red quasar fraction is consist ent with that required 
bv lGilli. Comastri. fc Hasingerl ()2007h to explain the quasar 
contribution to the X-Ray Background. 



4 TESTING THE KX METHOD 
4.1 Applying the KX method 

We tested the ability of our KX method variant to select 
quasars by constructing & b,j ~ R vs. R — K colour-colour 
plot in Figure [TTI In Figure [TT] the majority of the quasars 
are distinct from the stellar locus, only 7 lie within the stel- 
lar locus, successfully demonstrating the ability of the KX 
method to select quasars. At the top right of Figure 1111 the 
typical (mean) bj ~ R and R — K uncertainties due to pho- 
tometric uncertainties and quasar variability, are plotted as 
a cross. Examining the uncertainties in the quasar colours 
allowed us to determine how distinct from the stellar locus 
the quasars actually are. The photometric uncertainties in 
the bj a nd R magnitudes are both taken to be 0.1 mag- 
nitudes (jPrinkwater et al]|2000l ). The K magnitude photo- 
metric uncertainties were calculated by SExtractor during 
processing. To account for variability, we move both the bj 
and K photometry to the middle epoch of the r photome- 
try. To do this for the bj — R colour we used the structure 
function of iHook et al.l (| 19941 '). outlined in Section [3. II with 
Atobs = 15 years. We did not do this for the R — K colour 
as a structure function was not available for the K-band, 
instead the variability was estimated to be 0.1 magnitude s 
for all the quasars, from obs ervations in lEnva et al.l l|2002h . 
Data from lEnva et al.l (120021 ') was used as they measured the 
K band variability of quasars satisfying —29 < Mb < — 18 
and < 2 < 1, which closely matches the quasar sample 
constructed here. For the quasars in our sample with red- 
shifts greater than 1, 0.1 magnitude will be an over-estimate 
of the effect of variability; however, as it is of the order of 
the typical K magnitude photometric uncertainty, any over- 
estimate was negligible compared to the total uncertainty in 
the K magnitudes. Combining all the relevant uncertainties 
in quadrature the uncertainties in the bj — R and R — K 
quasar colours was calculated, and typical (mean) values 
plotted in Figure 1111 Examining Figure 1111 we concluded 
that the quasars in our sample were indeed distinct from 
the stellar locus; because, the bj — R and R — K uncertain- 
ties were small enough that the regions bounded by them 
were distinct from the stellar locus. 



4.2 Contrasting the KX method 

The argument for the use of the KX method in detecting 
quEisars is not that it finds red quasars; but that it is rel- 
atively bias free in its selection of quasars, whether red or 
blue. We experimentally tested this argument by compar- 
ing Figures [121 and [131 to Figure [141 Figure [121 is a standard 
U — bj vs. bj colour-magnitude plot used in selecting quasars 
in the UVX approach. Figure [131 is a [/ — 6j vs bj — R colour- 
colour plot used in the 2QZ to select quasars (|Croom et al.l 



BO 

a 

e 

Oi 

I 



— fT^ — — f^-- — V ' A ' — ^ 

\S \3 \4 \5 \6 \7 

\ V \ \ A A ^- 
\ \ \ \ -A \ 
\ A \ . A. A .. \ 
^ \ \ ■ v--':-.A:?f-^^^ ■ ■ \ - 
\ 1 \ \ \ "..A ... \ \ 

\ S^- ^\.::W;t"AA'XA' V 
\- ■ V ■ \ \ 

. \ \ ..Sv " 1\ \ o \ 

\ \ \-^.. \ \ \ - 

\ A iV ^ o\ \ \ 

\ \ A \ \ \- 

K ,\ \ \- 

\ A V V\ . \ „ \ 



\ \ \ \ 



_j , \ , , V. 



R-K (mag) 

Figure 11. bj — R vs. R — K colour-colour plot used to test the 
ability of the KX method to select quasars. Dots are all stellar 
morphology objects in the Fornax APMCAT catalogue. Solid tri- 
angles are blue {bj — k < 3.5) quasars and circles are red {bj — K ^ 
3.5) quasars in our sample. Where an R magnitude was unavail- 
able, limits for the bj — R and R — K colours are plotted using a 
square and arrow. Lines of constant bj — K are plotted as dashed 
lines, indexed as they increase from the bottom left corner. 
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Figure 12. U — bj vs bj colour-magnitude plot used to select 
quasars in the UVX method. Dots are all stellar morphology ob- 
jects in the Fornax APM catalogue with bj and U photometry. 
Solid triangles are blue {bj — K < 3.5) quasars and circles are red 
(bj — K ^ 3.5) quasars in our sample with U photometry. 
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Figure 13. U — bj vs b j — R colour-colour plot used to select 
quasars in the 2QZ survey llCroom et al Dots are all stellar 

morphology objects in the Fornax APM catalogue with bj and 
U photometry. Solid triangles are blue (fej — K < 3.5) quasars 
and circles are red (bj — K ^ 3.5) quasars in our sample. For 
the quasars in our sample without an R magnitude, squares arc 
plotted denoting a bj — R lower limit. 



Table 2. Example KX method criteria for selecting potential 
quasars for spectroscopic follow-up. 
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I2OO4I ). As Figures [121 and [131 only contain objects tiiat have 
U-band photometry, in Figure [T?] we generated the same 
colour-colour plot in Figure [TTI but only for objects that had 
U-band photometry so as to make a fair comparison. Com- 
paring Figures [T2] and [TSjto Figure [T^ qualitatively, the KX 
method immediately appeared more effective than the UVX 
method and comparable to the 2QZ multicolour method. We 
examined this quantitatively by measuring the percentage 
of quasars in our sample that would not be selected by the 
UVX and 2QZ multi-color methods, in total and as a func- 
tion of colour, and compared this to our KX method vari- 
ant. In Figure [T2I we use a colour criterion oi U — bj < -0.36 
to select quasars as an example of a t ypical UVX metho d 
criterion, based on Sra i th et al. (1997), iBovle et all l|l990l ): 
iLa Franca fc Cristianil (|l997l ). Using this cut 13 of the 69 
quasars, or 19%, were not selected in Figure [121 Breaking 
this into red and blue quasars, where red quasars satisfy 
bj — K ^ 3.5, 3 (10) of 11 (58) red (blue) quasars were 
not selected. The UVX method failed to select 27% of the 
red quasars compared to 17% of the blue quasars. In Fig- 
ure 1131 we applied the selection criterion used in the 2QZ 
l|Croom et al.l |2004| ). the criterion is marked by a dashed 



Figure 14. Repeat of Figure [TTI with the added criteria that only 
objects with U-band photometry arc plotted. 



line and objects below or to the left of the line were se- 
lected as quasars. This selection criterion failed to select 5 
of the 69 quasars, corresponding to 7% of the quasars. Bro- 
ken into red and blue quasars, 1 red and 4 blue quasars 
were missed, equating to 9% and 7%. For the KX method 
variant used here, in Figure [14] we applied an example of 
a suitable colour-based selection criterion, and plotted it as 
a dotted line in Figure [14] with objects below the line se- 
lected as quasars. The details of our example selection cri- 
terion are listed in Table [2] Our example selection criterion 
fails to select 7 quasars, 1 red and 6 blue, corresponding to 
10% of all the quasars and 9% and 10% of the red and blue 
quasars respectively. As we knew where the quasars were 
located on the colour-colour plot before defining a suitable 
colour-based selection criterion; comparing the results of our 
example KX method selection criterion to the results of the 
UVX and 2QZ methods is potentially an unfair comparison. 
As an alternative, we used the proximity to the stellar locus 
as a selection criterion i.e. whether a quasar was obscured 
by the stellar locus. In Figure 1141 7 quasars were clearly 
obscured by the stellar locus, 1 red and 6 blue. This cor- 
responds to 10% of all the quasars and 9% of both the red 
and blue quasars. Comparing all these percentages we found 
that the KX method variant was quantitatively superior to 
the UVX method, as the UVX method selected less quasars 
and was clearly biassed towards the selection of blue quasars. 
We also found that quantitatively the KX method variant is 
comparable to the 2QZ multicolour method in the number 
of quasars selected; however, the KX method variant is supe- 
rior as the 2QZ multicolour did demonstrate a bias towards 
the selection of blue quasars while the KX method variant 
did not. Our results are strong experimental support that 
compared to existing methods such as the UVX and 2QZ 
multicolor methods, the KX method, and in particular our 
variant, is ideal for selecting quasar candidates for spectro- 
scopic follow-up, independent of quasar colour. The key to 
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this result is that the lines of constant bj — K in Figure [14] 
cut through the stellar locus almost perpendicularly; there- 
fore, quasar bj — K colour does not affect the likelihood of 
quasar selection using our KX method variant. 



5 CONCLUSIONS 

In this paper, we constructed a sample of 69 quasars with 
measured bj — K colours, using a combination of APMCAT, 
IRIS2 imaging and FCSS spectroscopic identifications. 11 of 
these quasars are red, satisfying bj—K ^ 3.5, and all of these 
red quasars are radio-quiet according to the All Sky Opti- 
cal Ca talo gue of Radio/X-Ra y Sources. In accordance w ith 
iBarkhouse fc Halli (I2OO1I ') and lMaddox fc HewettI (I2006I '). as 
the 11 red quasars found here are unresolved, luminous and 
radio-quiet sources predominantly located at 2 > 1, this 
strongly indicates that the main cause of the red b,j — K 
colours is dust at the quasar redshift. Comparing our quasar 
sample to LBQS quasars, we found that the two bj — K 
colour distributions are different at the 99.9% confidence 
level, with our sample having a significantly broader, red- 
der distribution. Analysis of the uncertainties in the bj — K 
colours demonstrated that neither our red quasar detections 
or our comparison to the LBQS quasars is affected by them. 
A second analysis, on the effects of our datasets magnitude 
and completeness limits, revealed that they limited our mea- 
sured bj — K distribution to bj — K < 5. As red quasars are 
observed up to this colour limit, we concluded that the true 
bj — K distribution of quasars is even broader and redder 
than observed here. From the observed b,j — K colour dis- 
tribution, we robustly constrained the red quasar fraction 
of the K ^ 18.4 quasar population to be greater than 22%, 
and from a model estimated it to be 31%. 

Using the quasar sample constructed here, the viability 
of a KX method variant was tested. Using a bj — R vs R — K 
plot, the KX method variant was capable of separating the 
quasar sample out from the other objects in the APMCAT 
catalogue. Comparing the KX method variant to the UVX 
and 2QZ multicolour methods, the KX method variant was 
found to be as efi^ective in the number of quasars it selected; 
however, it was superior at selecting quasars independent of 
colour. For those reasons, the KX method variant used here 
is an ideal technique for future large surveys to use in select- 
ing potential quasars, whether red or blue, for spectroscopic 
follow-up. 
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